% 参数设置
am = 15;   % 调制信号的幅度
fm = 15;   % 调制信号的频率
fc = 50;   % 载波信号的频率
kp = 10;   % 相位调制指数

% 时间范围
t = 0:0.001:1;

% 调制信号
mt = am * cos(2*pi*fm*t);
subplot(3, 1, 1);
plot(t, mt);
title('调制信号');
xlabel('时间 (s)');
ylabel('幅度');

% 载波信号
ct = cos(2*pi*fc*t);
subplot(3, 1, 2);
plot(t, ct);
title('载波信号');
xlabel('时间 (s)');
ylabel('幅度');


% PM调制
pm_modulated = cos(2*pi*fc*t + kp*mt);
subplot(3, 1, 3);
plot(t, pm_modulated);
title('PM调制后的信号');
xlabel('时间 (s)');
ylabel('幅度');
